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PATENT APPLICATION 

TANGENTIAL VELOCITY MEASUREMENT USING 
INTERFEROMETRIC MTI RADAR 

STATEMENT OF GOVERNMENT INTEREST 

The United States Government has rights in this invention pursuant to Department 
of Energy Contract No. DE-AC04-94AL85000 with Sandia Corporation. 

FIELD OF THE INVENTION 

The present invention is generally related to radar systems. The present invention 
is also related to methods and systems used to determine a targets range and velocity. 
More particularly, the present invention is related to use of an interferometric moving 
target indicator radar to measure the tangential velocity component of a moving target. 
Multiple baselines, along with the conventional radial velocity measurement, allow 
estimating the true 3-D velocity vector of a target. 

BACKGROUND 

Although the present background describes the functionality and limitations of 
synthetic aperture radar systems or a particular class of communications, such description 
is merely provided to exemplify a problem capable of resolution with the present 
invention. Any discussion herein directed to specific radars or communications protocols 
should not be taken by those skilled in the art as a limitation on the applicability of the 
invention described herein. 

Radar systems use time delay measurements between a transmitted signal and its 
echo to calculate range to a target. Ranges that change with time cause a Doppler offset 
in phase and frequency of the echo. Consequently, the closing velocity between target 
and radar can be measured by measuring the Doppler offset of the echo. The closing 
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velocity is also known as radial velocity, or line-of-sight velocity. Doppler frequency is 
measured in a pulse-Doppler radar as a linear phase shift over a set of radar pulses during 
some Coherent Processing Interval (CPI). 

Radars that detect and measure target velocity are known as Moving-Target- 
Indicator (MTI) radars. MTI radars that are operated from aircraft are often described as 
Airborne-MTI (AMTI) radars. When AMTI radars are used to detect and measure 
ground-based moving-target vehicles, they are often described as Ground-MTI (GMTI) 
radars. 

In MTI radars, the angular direction of a target is presumed to be in the direction 
to which the antenna is pointed. Consequently, a MTI radar generally offers fairly 
complete position information (angular direction and range) to some degree of precision, 
but incomplete velocity information since Doppler is proportional to the time-rate-of- 
change of range, i.e. radial velocity. Tangential velocities, that is, velocities normal to 
the range direction do not cause a Doppler shift, so are not measured directly. Tangential 
velocities can be measured indirectly by tracking the angular position change with time, 
but this requires a somewhat extended viewing time for any degree of accuracy and/or 
precision. 

Multiple MTI systems might be employed in concert, each measuring radial 
velocities in different spatial directions. In this manner, a two-dimensional (or even full 
three-dimensional) target velocity vector may be estimated. This technique, however, 
requires that the radars be widely separated to facilitate the necessary triangulation (e.g., 
being based on different aircraft in the case of GMTI systems). 

GMTI systems are often employed from moving radar platforms such as aircraft, 
that is, the radar itself is in motion with respect to the ground. Consequently, the 
stationary ground itself offers Doppler frequency shifts. In addition, since different areas 
of the ground are within view of different parts of the antenna beam, and have somewhat 
different radial velocities, the ground offers a spectrum of Doppler frequencies to the 
radar. This is often referred to as the clutter spectrum, and can mask the Doppler returns 
for slow-moving target vehicles of interest. Of course, if a target's Doppler is outside of 
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the clutter spectrum, its detection and measurement are relatively easy. This is called 
"exoclutter" GMTI operation. Detecting and measuring echo responses from slow- 
moving target vehicles that are masked by the clutter are considerably more difficult, and 
are called "endoclutter" GMTI operation. 

The ability to observe targets masked by clutter is often called "sub-clutter 
visibility." Reducing the effects of clutter on detecting and measuring such targets' 
motion is often termed "clutter suppression." This is most often accomplished by 
employing multiple antennas on a single aircraft arrayed along the flight direction of the 
radar, and is often called a Displaced Phase Center Antenna (DPCA) technique, or 
Interferometric GMTI. The following patents provide background information on the use 
of more than one antenna in radar systems: US Patent 4,885,590, issued December 5, 
1989 to M. A. Hasan, entitled "Blind speed elimination for dual displaced phase center 
antenna radar processor mounted on a moving platform"; US Patent 4,086,590, issued 
April 25, 1978 to W. B. Goggins, entitled "Method and apparatus for improving the 
slowly moving target detection capability of an AMTI synthetic aperture radar"; US 
Patent 5,559,516, issued September 24, 1996 to J. A. Didomizio, R, A. Guarino, entitled 
"Dual cancellation interferometric AMTI radar"; US Patent 5,559,518, issued September 
24, 1996 to J. A. DiDomizio, entitled "Low target velocity interferometric AMTI radar"; 
and US Patent 5,818,383, issued October 6, 1998 to E. F. Stockburger, H. D. Holt Jr., D. 
N. Held, R. A. Guarino, entitled "Interferometric moving vehicle imaging apparatus and 
method." 

Interferometric techniques allow making independent angle measurements not 
affected by target motion, thereby facilitating discrimination of a moving vehicle in one 
part of the antenna beam from clutter in another part of the antenna beam that otherwise 
exhibits identical Doppler signatures. Interferometers can be constructed from separate 
distinct antennas, or from monopulse antennas that offer the equivalent of separate 
distinct antenna phase centers in a single structure. Although interferometric systems and 
method provide improved target analysis through clutter reduction, the acquisition of 
tangential velocity measurements within a single CPI still remains problematic, and has 
not been adequately addressed in the art. There remains a need for a more complete 
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target velocity vector measurements and analysis for time-critical moving vehicles. 
Furthermore, there is a need that such measurement and analysis continue to be provided 
from the vantage point of a single system. 
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SUMMARY OF THE INVENTION 

A radar interferometer can measure angular position to a target with a great deal 
of precision, even with a single radar pulse. It does so by measuring the phase difference 
between echoes arriving at the two antennas. A target with tangential velocity will 
exhibit a pulse-to-pulse change in the angular position as measured by the phase 
difference between the antennas. This manifests itself as an interferometric phase that 
changes with time, i.e., a Doppler difference frequency. By measuring this Doppler 
difference frequency over some Coherent Processing Interval (CPI), a tangential velocity 
can be calculated for the target. This tangential velocity will be in the direction of the 
interferometric baseline. Consequently, a multiple orthogonal baseline arrangement can 
measure tangential velocities in both the azimuth and elevation directions. These coupled 
with the radial velocity derived from traditional Doppler processing enables a full 3- 
dimensional velocity vector to be measured from a single CPI. 

According to features of the present invention, subaperture techniques allow for 
filtering individual Doppler returns when multiple moving targets exist at the same range. 

The techniques of the present invention are usable for a wide variety of radar 
systems applications, including air traffic control, ground vehicle target tracking, law- 
enforcement, and traffic monitoring and control. Techniques of the present invention 
also extend to other coherent remote sensing systems such as sonar, ultrasound, and laser 
systems. 



6 



Attorney Docket No. 1000-1298 



PATENT APPLICATION 

BRIEF DESCRIPTION OF THE DRAWINGS 

The accompanying drawings, which are incorporated in and form part of the 
specification, illustrate embodiments of the present invention and, together with the 
description, serve to explain relevant principles of the invention: 

Figure 1 provides an illustration of a vector r, from target exhibiting radar echo to 
radar, where r c is a vector from target reference location to radar, s is a vector from 
target reference location to target exhibiting radar echo, andr, = r c -s ; 

Figure 2 provides an illustration of a geometry where a second receive-only 
antenna is located with a baseline vector b, with respect to the original transmit/receive 
antenna, and the bistatic return echo exhibits a video phase; 

Figure 3 illustrates method steps for processing Exoclutter GMTI signals in 
accordance with features of the present invention; 

Figure 4 illustrates additional method steps for processing Endoclutter GMTI 
signals in accordance with features of the present invention; 

Figure 5 illustrates method steps for operating a stationary MTI system model in 
accordance with features of the present invention; 

Figure 6 illustrates method steps for operating an airborne GMTI system model 
in accordance with features of the present invention; and 

Figure 7 illustrates a block diagram of a system and its components in accordance 
with features of the present invention. 
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DETAILED DECRIPTION OF THE EMBODIMENTS 

Although the following provides a detailed description that includes examples of 
the functionality and limitations that can be found in radar systems or may refer to a 
particular class of signal detection and assessments using radar, such examples are 
merely provided to exemplify a signaling problem commonly found in the radar art. Any 
discussion herein directed to specific radars and their associated protocols should not be 
taken by those skilled in the art as a limitation on the applicability of the invention 
described herein. 

With a monostatic radar employing a Linear-FM chirp, stretch-processing, and 
quadrature demodulation, it is well-known that the return echo exhibits a sampled video 
phase that can be adequately described by: 

®v{Un)=l(<v + ?T s i%r c \-\r s \) (1) 
where, 

c = velocity of propagation, 
co = radar center frequency, 
y= radar chirp rate, 

T s = Analog to Digital conversion period (sampling period), 

i = intra-pulse sampling index with - 1/2 < i < 1/2 , 

n = inter-pulse index with - N/2 <n< N/2 , 

r c = vector from target reference location to radar, 

r s = r c ~ s = vector from target exhibiting radar echo to radar, and 
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s = vector from target reference location to target exhibiting radar echo. 

The geometry for this expression is illustrated in FIG. 1 wherein a monostatic radar 101 
and target 102 are identified. As shown in FIG. 1, vectors r c and s can be presumed to 
be able to change with index n. The signal itself will have some amplitude A, and with 
this phase can be described by: 

X v (i,n)= Ae j<S> ^' n \ (2) 
The phase is adequately approximated by: 

^v(hn)=% + 7 T s i) r "°" 



r c\ 



(3) 



Referring to the geometry illustrated in FIG. 2, if a second receive-only antenna 211 is 
located with a baseline vector b with respect to the original transmit/receive antenna 101, 
then its bistatic return echo from target 102 exhibits a video phase described by: 



®v,b(i,n) = -(a>+yT s i%r c +(r c + bj-|r, + {r s +bj), (4) 
which for small baseline lengths can be adequately approximated by: 



2/ „ .\( r, 0 s bos 



c I l r d 2 Yc\j 



(5) 



(6) 



The video signal that exhibits this phase can be described by: 

X v ^(Un)=Ae^ n \ 

When |b| = 0 , the monostatic case is reduced. 

A moving target implies that target location vector s changes from radar pulse to pulse as: 
s = s 0 +\ s T p n (7) 
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where, 



s 0 = target reference position at n = 0, 
v s = target velocity vector, and 
T p = Pulse Repetition Interval (PRI). 
In general the radar itself can be moving, i.e. changing position from pulse to pulse as: 



r c=r c0 +y c T p n 



(8) 



where, 



r 0 = radar reference position at n = 0, and 
v c = radar velocity vector. 

Assuming a constant baseline vector, an expansion to the approximation can be provided 
as: 



(( 



r c0° S 0 b °Sn 
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(9) 



Typical GMTI parameters are such that this can usually be further simplified to the an 
approximation, as follows: 



2/T s 



V l r co| 



b°s 0 

2|r 



coiy 



i + 
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+ ■ 
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\ 1 
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c0\ J 



(10) 



Some constant phase terms that are inconsequential to target position and motion 
estimation have been intentionally ignored with the foregoing. 
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Given the foregoing analysis, the phase is then of a form: 
^K.b ('>")= <*>ii + o) n n + <p 
where, 



(11) 



o>i = 



<»n = 



2yT s 



r c0° s 0 | b °Sp 

< | r co| 2 \r c o\j 



2coT„ 



V l r co| | r co' 



2\r, 



cO 



, and 



2a> 

<P= — 
c 



There are three principal terms for the present analysis. The first is a coefficient 
of index i, and represents a frequency with respect to index i. The second is a coefficient 
of index n, and represents a frequency with respect to index n. The third is a constant 
phase term with respect to indices i and n. A complete MTI data set corresponds to a 
data array of size Nx I, with elemental positions denoted by indices n and /. 

The coefficient of index i can be identified via a Fourier Transform of the video 
signal data over index i. The coefficient of index n can be identified via a Fourier 
Transform of the video signal data over index n. Performing Fourier Transforms across 
both indices i and n yields a range-Doppler map, where peak responses in the map 
indicate targets at a specific range and velocity. That is, the range-Doppler map can be 
described by: 

Z b {v,u)= 3„(3 i (^(/,«))) = YL X v( i > n Y im '* Jm '* = A ^vk - 



(12) 



where, 



v = the range index of the range-Doppler map with - V/2 < v < V/2 , 
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u = the Doppler index of the range-Doppler map with -U/2 < v < U/2 , 

o v = the map's range scale value corresponding to index v, 

co u = the map's Doppler scale value corresponding to index u, 

W v { ) = the map's range impulse response function, with peak at W v (o), and 

w u ( ) = the map's Doppler impulse response function, with peak at W u (o) . 

The range-Doppler map for the monostatic case can now be defined as: 

z oM=z b (v, w )| b=0 . (13) 

A coordinate frame with unit vectors can now defined as: 
f = - r co/| r co| = ran ge direction, such that r c0 = -|r c0 |? , 

a = azimuth direction vector, horizontal, to the right, and orthogonal to f , and 
e = a x f = elevation direction vector. 

Tangential velocities are in the azimuth and/or elevation directions and are 
embodied in the b o \ s term. 

The vector quantities in this frame can also be defined as: 

r co =-Ko|r 
v c =v c ,f + v ca a + v ce e 
s 0 =s r r + s a a + s e e 

b = b r r + b 0 a+b e e (l4) 
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A more simplified expression is therefore possible as follows: 



bos. 



2r 



2a)T 
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(15) 



Exoclutter GMTI 

Exoclutter GMTI generally is monostatic (i.e., no baseline is relevant). Furthermore, it 
assumes that (or resigns itself to the case where) target radial velocity offers Doppler 
shifts much greater than the clutter spectrum width, namely |v Jr | » |v c os 0 |/|r c0 |. 
Consequently the model for target phase becomes: 

*F,bM= — (-*r> + ^-(-V Jr )«. (16) 

c c 

Referring to FIG. 3, a flow diagram providing steps for processing the Exoclutter 
GMTI signal is illustrated. According to the method, CPI is collected 310. A range 
transform is then performed 320. After the range transform is performed, the system can 
Identify a target's range position 330. Regardless if range position is first performed as 
shown in step 330, the system also performs a Doppler transform 340. After Doppler 
transform is conducted 340, a target's radial velocity can be determined. 350. 

Target radial position is measured by a range transform across index z, and target radial 
velocity is measured with a Doppler transform across index «. The range-Doppler map 
becomes approximately: 



Z 0 M= A W v [co v -^(-* r )VL ~^(-v jr ) 
\ c ) y C ) 



(17) 
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Endoclutter GMTI 

Endoclutter GMTI uses interferometry with a baseline, and allows for measuring target 
radial velocities with Doppler shifts less than the clutter spectrum width, namely 
Kr|<| v c ° s o|/| r co|- ^e baseline can be assumed to be small enough that it doesn't 
influence significantly the result of the Doppler transform across index n. Furthermore, 
the baseline should be generally aligned in the azimuth direction, and horizontal radar 
flight path is presumed. Consequently, the model for target phase can be: 



which is still in a form of equation (11), namely O k b (/, n) = co .i + <o n n + <p . 

Referring to FIG. 4, the processing steps for the Endoclutter GMTI signal are shown. As 
shown in Step 405, CPI is collected from a first antenna, while CPI can also be collected 
using a second antenna 407. Next, a range transform can be performed 410 and 413 
using information from both antennas. A targets range position can be identified 415 at 
this point. Next Doppler transform is concurrently performed 420, 423 by both antennas. 
Range correction, if applicable, can then be performed at both antennas at a radar as 
shown in step 430, 433. The system can then perform interference 440 given data 
provided from both antennas, after which target azimuth positions can be identified 445. 
If Target azimuth position requires correction, correction can occur as shown in step 450. 
Once the range and azimuth positions are determined, a radar system in accordance with 
the present invention can identify a target radial velocity as shown in step 460. 

A range transform across index i enables the identification of target radial position 
s r . That is, 




(19) 



where the remaining phase is: 
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®R,b( v >")= °> n n + <P- ( 20) 

A Doppler transform across index n identifies <y„, which in turn identifies the quantity 
( v sr ~ v cr ^/|r c0 | - v ca 5 a /|r c0 |). That is, the range-Doppler map is described by: 

^M=^ n (x Rb {v,n))=A W v {fo v -(o^W u {co u -co n yP . {2 \) 

Since s r is known, its effects can be compensated to identify more specifically 
( v sr-v ca s a /\r c0 \), but still leaves some ambiguity in resolving energy to some 
combination of the unknown quantities radial target velocity v sr and target azimuth 
position s a . The monostatic antenna does not contain the baseline term in the above 
model, but the bistatic second antenna does contain this term. 

The range-Doppler maps made from both ends of the baseline can now be 
compared. The two results are compared by interfering them to yield: 

IF(v,i)=Z b (v,u)Z 0 *(v,u) = \A W v {o> v -c>„f e J ^ v > u \ (22 ) 
where * denotes complex conjugate, and at each pixel: 



A**(v,h) = — 
c 



\ 2 \ r cO\J 



(23) 



Identifying this phase enables an independent estimate of target azimuth position s a , 
which can then be used to correct the result of the earlier Doppler transform result to 
uniquely identify target radial velocity v sr . These quantities can thereby been separated, 



or discriminated. 



It is important to understand that the baseline's influence is analyzed only after the 
complete range-Doppler maps are formed for both receiving antennas. 
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Stationary MTI for Estim ating Tangential Velocity 



While extracting tangential velocity information, use of a stationary interferometric radar 
(v c = 0) with baseline oriented in the azimuth direction can be considered. The model for 
target phase determination is: 



-s.+- 



b.s. 



2r 



cO\J 



2coT 



-v.. + 



b v 

a r sa 



n + 



c0 J J 



2a 
c 



b^ 
2lr 



(24) 



Referring to FIG. 5, the processing steps for operating a stationary MTI model are 
illustrated. According to the method, CPI is collected from first and second antennas as 
shown in blocks 510 and 520. After CPI collection, range transform is performed 515, 
525, where after a target's range position can be identified 527 given range 
transformation data from the first antenna; and additionally interference can be performed 
529, given range transform data from both antennas. Range transform data from the first 
antenna can be used to perform Doppler transform as shown in block 530, where after a 
target's radial velocity can be identified 535. After Interference is performed 529, a 
targets azimuth position can be identified as shown in block 540. A target's azimuth 
velocity can also be identified as shown in block 555, after Doppler transform is 
performed in block 550. 



A range transform across index i enables the identification of a target's radial 
position s r , and also leaves the data with a residual phase: 



®R,b = 



2<oT/ 



~ v sr + 



b a Vsa 



2r, 



c0| ) 



2a> 
c 



b a s a 



2r, 



c0 



(25) 



Target radial velocity can be estimated from the monostatic antenna data in the usual 
manner. If, however, this range-compressed data from each of the two antennas is 
interfered with each other, that is, a phase comparison is made before any Doppler 
transform, the resultant data can be expressed as: 



IF R {v,n) = \AW(o) v -o n fe j ^ v ' n \ 



(26) 
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where, 



A*(v,n) = p - 



f 7 \ 

b a v sa 

V 2 l r co| y 



n + 



2o)\ b Q s a 



(27) 



An interference signal with this phase characteristic is generated by point-by-point 
multiplication of the data from one range-compressed data set with the complex 
conjugate of the data from the other range-compressed data set. 

The coefficient of index n in the phase is now a Doppler difference frequency that 
depends on target azimuth velocity v SQ . That is, for the interference signal now: 



2o)T n 



co n =• 



^ 2 | r co|; 



(28) 



A Doppler transform of this interference signal over index n can now enable 
identification of target azimuth velocity v sa corresponding to the frequency content of the 
interference signal; that is, the range-Doppler map for this interference signal is now 
described by: 



Zjf (v, u) = 3„ (IF R (v, n)) = \A W v (a> v -co t f W u (co u - co n )e^^ , 



(29) 



whereby the 2-dimensional peak can now describe target range and target tangential 
velocity. The average phase of the interference signal remains dependent on target 
azimuth position s ay which is: 



A*(v,0) = ^ 
c 



b a s a 
2 \ V co\) 



(30) 



Regardless, the target tangential velocity in the azimuth direction v sa can also be 
identified. 



The shortfall of the stationary MTI model is that it can be very sensitive to noise, because 
data that is interfered is only range compressed up to this point, and doesn't benefit from 



17 



Attorney Docket No. 1000-1298 



PATENT APPLICATION 

the noise reduction offered by Doppler processing. Furthermore, multiple targets at the 
same range but at different radial velocities are indistinguishable from each other, and 
can in fact severely diminish the ability of a system to find the correct tangential velocity 
for any target. Nevertheless, the concept of tangential velocity derived from 
interferometric MTI can now herewith be established. 



More Robust Airborne GMTI Radar for Estimating Tangential Velocity 

A more complex scenario involving a moving radar can be addressed using 
airborne GMTI radars. For simplification of the following example, elevation velocities 
have been limited to zero, and baseline orientation to the range-azimuth plane. 
Furthermore targets of interest will be assumed to exist in the exoclutter region. 

The complete phase model for the data has been determined as: 



2yT s 
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(31) 



Referring to FIG. 6, the processing steps for operating a airborne GMTI radar are 
illustrated. According to the method, CPI is collected from first and second antennas as 
shown in blocks 610 and 620. After CPI collection, range transform is performed 615, 
625, where after a target's range position can be identified 617 given range transform 
data from the first antenna. Range position can be corrected as shown in blocks 630 and 
640, where after Doppler transform can be performed within subapertures as shown in 
blocks 650 and 660. A target's radial velocity can be identified 655 given the Doppler 
transform data from block 650. Radial velocity of the target can be corrected as shown in 
blocks 670 and 680, where after interference can be performed 685. The targets azimuth 
position can be determined 687 after Interference is performed 685. Furthermore, any 
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phase shift between subaperatures can be identified as shown in block 695, from which a 
target's azimuth velocity can be identified as shown in block 695. 

A range transform of the data across index i enables the identification of a target's 
radial position s r , and leaves the data with a residual phase: 



v s +v s by +6 v 

y + cr r ca^a r sr Ly a y sa 



"cO 



2 k 



2co 



c0\ J 



b r s r +b a s a 



2r, 



cO 



(32) 



A function exhibiting some phase 0 perturbed by an undesired but known phase s can be 
corrected by multiplying with a phase correction signal of unit amplitude and the 
negative of the phase perturbation. Accordingly: 



(Ae j(Q+e) )e- Je = Ae i{@+eyi£ = Ae i@ . 



(33) 



In this manner, the data can be corrected for its influence by applying a phase correction 
to yield: 



2a 



®R,b( v > n ) = ®R,b( v > n ) 



^ b r s r ^ 



^ 2 i r c0|y 



(34) 



because target radial position s r is now known, or more explicitly, 
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(35) 



At this point the CPI is split into two subapertures by dividing along index n to yield two 
new indices m and k such that: 



n = m + \ k — , 

I 2)2 



(36) 



where within a subaperture - N/4 <m< N/4 , subaperture index k takes on values 0 or 1 . 
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The range compressed data can then be modeled with exhibiting phase as shown in steps 
630 and 640 of FIG. 6, according to the following: 



®'R,b( V > m > k ) = 
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(37) 



Given exoclutter targets as the specimen, a Doppler transform across index m yield: 
estimate of target radial velocity v sr and a residual phase term: 
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(38) 



The data can be corrected for its influence because a target radial velocity v, r is now 
known by applying a phase correction to yield: 



®'D,b( V > U > k ) = 



(oT p N 



2a) 



bgS a 

\A r co\) 



\ \ r c0 
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2 k, 



c oi;j 



(39) 



The partially processed data from the two antennas can then be interfered with on a point- 
by-point basis. This would yields an interference signal result with phase given by: 



aT n N(bv A 
u a v sa 



, 2a 
k + 



b a s a ) coT n N 



y 2 \ r co\ 



2c 



bgV S g 

k 2 Ko|. 



(40) 
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PATENT APPLICATION 

This interference result still has two subapertures. The average phase allows estimation 
of target azimuth position s a . Compare the phase of corresponding interference data 
across the two subapertures, e.g., across index k, a phase difference measure of: 



A(A^(v,«,/:)) = A^(v, W ,l)-A^(v, W ,0)= 



b a v sa 



2 k 



[ c0 



(41) 



can result, which allows unique identification of target azimuth velocity v sa ; which 
means, the phase difference at each pixel of the corresponding subaperture interferograms 
is proportional to tangential velocity v sa . 

In principle, splitting the CPIs into more than two subapertures prior to interfering them 
also allows extraction of tangential velocities. Furthermore, it can now likely be 
appreciated that a third receive-only antenna located with a second baseline vector 
oriented in the elevation direction would allow the additional discerning of the elevation- 
direction tangential velocity in a similar manner. 



Accuracy and Precision of a Measurement 

Interferometric measurements tend to be very sensitive to noise, and work best for targets 
with relatively high Signal-to-Noise Ratio (SNR). In this case, a phase difference can be 
measured to within an RMS error given by: 

aM= im- < 42 > 

and a radian frequency can be measured to within an RMS error given by: 

a "'Tm- (43) 

where Tis the RMS observation interval. 
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PATENT APPLICATION 

Consequently, for the more robust airborne GMTI as described above, including two 
subapertures per antenna, the precision of the tangential velocity measurement can be 



shown to be: 



CT v„ =1- 



[ c0 



\ b a J 



where SNR refers to the results of the interference operation. 

Referring to FIG. 7, a callout from an airborne MTI radar system 700 that is 
illustrated in operation tracking a ground-based vehicle 705. According to the call-out, 
the radar system 700 includes a first antenna 710 and second antenna 720. The antennas 
provide acquired signals to the modules comprising the MTI. Modules can include 
Range Transform 730, Doppler Transform 740, Interference 750 and Phase shift identity 
760 capabilities. An airborne radar system with modules as depicted in FIG. 7 can 
support the methods described herein and illustrated in FIGS. 3-6. 



Example: 

An interferometric GMTI has an antenna baseline spacing b a = 3 m, and an operating 
frequency of 16.7 GHz. Furthermore, a target exists with an actual tangential velocity of 
90 m/s at a range |r c0 | = 5 km, and SNR after interference is 100 (20 dB). With a CPI of 
TpN= 0.25 s, an expected system precision of = 3.8 m/s can be expected. 

The result of the previous example assumes a single CPI. It should be appreciated after 
the foregoing teaching that tracking the target over multiple CPIs in the conventional 
manner will allow refinement. 

The embodiments and examples set forth herein are presented in order to best 
explain the present invention and its practical application and to thereby enable those 
skilled in the art to make and utilize the invention. However, those skilled in the art will 
recognize that the foregoing description and examples have been presented for the 
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PATENT APPLICATION 

purpose of illustration and example only. The description as set forth is not intended to 
be exhaustive or to limit the invention to the precise form disclosed. Many mogadons 
and variations are possible in light of the above teaching without departing from the spirit 
and scope of the following claims. 
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